On the Use of Finite-Size Scaling to Measure Spin-Glass Exponents 
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Finite-size scaling (FSS) is a standard technique for measuring scaling exponents in spin glasses. 
Here we present a critique of this approach, emphasizing the need for all length scales to be large 
compared to microscopic scales. In particular we show that the replacement, in FSS analyses, of 
the correlation length by its asymptotic scaling form can lead to apparently good scaling collapses 
with the wrong values of the scaling exponents. 



This paper is a critique of finite-size scaling (FSS) 
methods as applied to spin glasses. We focus primar- 
ily on one- and two-dimensional Ising spin glasses, which 
order only at temperature T — 0, though the underly- 
ing ideas are quite general. The model is defined by the 
Hamiltonian H = — X]<i j> JijSiSj — hj^i^i, where h 
is an external field which will be zero unless stated oth- 
erwise. 

In dimension c? = 2, a range of different methods for 
measuring exponents has led to widely differing quoted 
values, sometimes by as much as a factor two, for the 
same exponent Qllllllli, HHIl. Despite years 
of numerical study, these discrepancies have never been 
satisfactorily resolved. Here we present analytical and 
numerical results in dimension d = 1, for which exact 
analytical values for the scaling exponents are known. 
By applying standard FSS methods to the data, we find 
large discrepancies between the numerically determined 
exponents and the exact ones. We show that these dis- 
crepancies disappear if the data are plotted in a different 
way, using the exact correlation length in the FSS analy- 
sis instead of its leading scaling form. We conclude that 
the main cause of the discrepancies is large corrections to 
scaling in the expression for the correlation length, rather 
than corrections to FSS itself. 

For the critical behavior as T — s- 0, there is only one in- 
dependent exponent JS| , except perhaps when the ground 
state has a non-trivial degeneracy. The exponent on 
which we focus our attention is the 'stiffness exponent' 
6, which describes the dependence on length scale, I, of 
the energy, E, of an excitation from the ground state: 

~ flfiESini- For d < 2, 6* is negative and large 
excitations are easily created by thermal fluctuations, de- 
stroying the ground-state order. Setting E ^ T gives 
a characteristic length scale ^ ~ T~'^ , with v — —1/9, 
above which the ground state is unstable against thermal 
fluctuations, i.e. ^ can be identifled with the correlation 
length [l[l3. 

The exponent describes the dependence on I of the 
energy of a 'droplet' of reversed spins, of linear size I. The 
energy is associated with the boundary of the droplet. A 
measurement of 9 can be made through the numerical 
study of domain wall energies, in which a wall, or 'in- 
terface', is imposed on a finite-size system by choice of 



boundary conditions 1]. One finds 9 k, —0.28 in d = 2. 
It was recently shown that this result can be obtained 
in a way which is explictly independent of the particular 
boundary conditions used to impose the interface 12!|. 
The scaling result v = —1/9 then predicts w 3.6 for 
the correlation length exponent. 

There has been some debate as to whether the expo- 
nent 9 (sometimes called 9dw) extracted from domain 
wall energies is identical to the exponent obtained by 
studying directly the size dependence of droplet energies. 
Droplet excitations can be created from the ground state 
by reversing a central spin holding the boundary spins 
fixed, where the scale / is the linear size of the system 
[l3| . This approach gives 6 « —0.42, though the analysis 
has been criticised by Middleton who argues that only 
droplets whose area exceeds some fraction of P should 
be included Hartmann and Moore |0| have shown 
how the apparent difference between the two values of 
9 can be explained by invoking a correction to scaling 
(a subdominant term of order , with uj > —9, in the 
expression for the droplet energy), and that the corre- 
sponding value of 9 is consistent with that obtained from 
domain wall studies. Such corrections to scaling may also 
account for the value 9 ~ —0.46 obtained by Picco et al. 
in their recent study of droplet energies Q . 

A 'direct' measurement of the exponent v, that de- 
scribes the divergence of the correlation length, can be 
made by studying, for example, the spin-glass suscepti- 
bility 



XSG 



(1) 



where the angular brackets represent a thermal average, 
the overbar is a disorder average, and L is the linear size 
of the system. At T = every term in the sum is unity 
and XSG = L^- For T > 0, Js's])^ = fir/0, where 
f{x) is a scaling function with /(O) = 1, f{x) ~ exp(— a;) 
for large x, giving '9] XSG ~ ~ T^^ with 7 — dv. 
FSS predicts xsG = L'^F{L/S,) for L 00, ^ ^ 00 
with L/i fixed but arbitrary. Using ^ - T-" = T^/\ 
one can determine 9 by plotting L~'^XSG against TL^^, 
and choosing 9 to give the best data collapse. Using 
this method, Kawashima et al. |^ find i' « 2.0, i.e. 9 « 
—0.5, which differs significantly from the value ~ —0.28 
inferred from the domain- wall studies. 
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A second method that has been used to extract 9 
numerically is to measure the magnetization per spin, 
m{h) = L~'^_y2ASj), induced by a small magnetic field 
/i at T = Id lia 111 III 113 Since in zero field we have 
m(0) ~ L~'^/^, a simple scaling argument for small h 
gives m{h) = L-'^/^g[h[L)/J{L)], where h{L) ~ hL'^''^ is 
the effective field at scale L and J{L) ~ is the effective 
coupling at this scale. Hence 'm{h) = L-'^^^M{L'^/^-'^h). 
In the thermodynamic limit m should be become inde- 
pendent of L, which implies m ~ h^^^ with d = l—20/d w 
1.28 for d = 2. An alternative way of writing the FSS 
form is 



1.5 



rriL 



(2) 



with G(0+) = const and G(x) ~ x for x — > oo. With 
this method, Rieger et al. |a] obtained S w 1.48, which 
differs significantly from the scaling prediction « 1.28 
and, naively, predicts that 6 w —0.48 instead of the value 
« —0.28 obtained from domain-wall studies. 

To summarise, different ways of measuring 9 give dif- 
ferent, and seemindy incompatible, results. The value 
1/ « 2 obtained in j5j is equivalent to 6* w —0.5, which 
is similar to the value 9 « —0.48 inferred from the m{h) 
data of and to the values w —0.42 and « —0.46 ob- 
tained from studying droplet excitation in and 0] 
respectively, all of which differ from the value « —0.28 
obtained from domain- wall studies 0, 0, ^3 ■ 

In an attempt to understand these differences, we have 
carried out analytical and numerical studies in space di- 
mension d = 1, for which the corresponding exponent 
values are known exactly. We mimic the two-dimensional 
studies of Kawashima et al and Rieger at al 0], and 
look at the T-dependence of XSG at h = and the h- 
dependence of m at T = respectively. We conclude 
that these that these quantities are affected by correc- 
tions to scaling so large that it is essentially impossible 
to extract the correct exponent values from system sizes 
that are accessible in d = 2. 

In the remainder of the paper, therefore, we consider 
the d ~ 1 Ising spin glass in two situations: (i) T > 
and h ~ 0. We calculate xsg and use the FSS form 
XSG — Lf{TL~^) to determine 9. We show that the 
exact value of 9 gives a very poor data collapse for the 
system sizes studied, and that a reasonable data collapse 
is obtained with a significantly different value of 9; (ii) 
T ^ 0, h > 0. We calculate and use the FSS 

form mL{h) ~ L^/'^g{L^/'^h^/^) to determine 5. Again, 
the exact 6 gives a poor collapse, and the best collapse is 
obtained with a very different 5. To facilitate comparison 
with the d = 2 data of Kawashima et al and of Rieger et 
al. we choose, for the h = results, a bond distribution 
P(J) engineered to give 9 — —0.282, i.e. a value equal 
to that of the d — 2 system, while for the h > data 
we choose the distribution such that 6 = 1.282 is equal 
to the value predicted in d = 2 for 9 — —0.282. These 




FIG. 1: Susceptibility scaling plot. From top to bottom: 
—0.28 (theoretical) and 6 = —0.36 (best fit by eye). 



choices are are imposed by using a bond distribution of 
the form 



Pa(J) oc |J|"exp(-jV2) 



(3) 



For any distribution satisfying P{J) ^ | J|" for J ^ it 
may be shown 1^ that 9 = -1/(1 -I- a) for 0? = 1, so the 
choice a = 2.546 gives 9 = —0.282, while a = 6.042 gives 
S = 1.282. 

(i) T > 0, = 

The Ising Hamiltonian for li = 1 can be written H = 
— J2i JiSiSi+i, and we use free boundary conditions. It 
is straightforward to show that (for j > i) (SiSj) = 
n^~}tanh(/3J,) and 



(4) 
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FIG. 2: Data collapse with scaling variable L/£^exac 



where a = tanh {(3Jr) is independent of r. Finally 



XSG 



l + a 2a(l-a^) 



1 



L{1 



(5) 



The quantity a — J dJ Pq, (J) tanh^(/3 J) may be eval- 
uated numerically for any temperature T. The result- 
ing XSG is plotted in Fig. 1, in the scahng form xso/L 
against TL~^, for lattice sizes up to L = 320. In the 
upper figure, we use the exact value, 6 = —0.282, while 
in the lower we use 9 = —0.36, which is our 'best fit by 
eye'. Comparing these two plots one observes that (i) 
The exact value of 9 gives a good collapse only at small 
values of the scaling variable, where xsg/L is close to 
its T = value of unity; (ii) The 'best fit by eye' is a 
much better fit over a large part of the plot, especially 
for the larger systems. Small but systematic departures 
from perfect scaling are evident in the low T region, but 
these are only observable because we have perfect data 
(no statistical errors). In real (i.e. noisy) data such small 
effects could easily be obscured by the noise, and might 
lead one to suppose that the correct value of 9 were close 
to -0.36. 

This tell us that the relation ^ ^ T^"^, where v = 
— 1/9, is only valid for rather small T. Corrections to 
this form are important over most of the regime presented 
in Figure 1. A striking confirmation of this is provided 
by Figure 2, where the same data are plotted against 

exact 1 

where iexact is the exact correlation length for 
each temperature. From Eq. (2) we can identify this 
length scale as S,exact = — 1/lna. The data in Fig. 2 
collapse almost perfectly for all sizes L > 40. 

The lesson here is that it is not FSS itself which is 
breaking down, but the use of the relation ^ T~'^ over 
the whole temperature range explored. We suspect that 




ln[h 7mL(h=0)] 



FIG. 3: Magnetisation versus field, scaled by the numerically 
determined zero-field magnetization. Top: theoretical, 5 — 
1.282 and Bottom: best collapse by eye, Scoi = 1.37 



similar problems affect the interpretation of the data of 
Kawashima et al. 01 . 

(ii) r = 0, h > For non-zero magnetic field, the 1-d 
problem cannot be solved analytically in closed form for 
general system size, L. One can, however, determine the 
relevant correlation length, ^, numerically. Here ^ is to 
be interpreted as the length scale over which the ground 
state for h > locally resembles the h ~ ground state. 
We define it as the average 'domain length', where a do- 
main is a cluster of spins completely aligned with one 
of the two h = ground states. This definition corre- 
sponds to the scaling argument L'^/'^h^/^ in Eq. (|2Jl, i.e. 
this argument - (L/^y^^ with ^ oc ft-^^/d^^ 

Data for system sizes up to L = 680 are displayed 
in Fig. 3, where each data point represents an average of 
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FIG. 4: Magnetization scaling plot using the numerically de- 
termined correlation length 



50 000 samples. As explained above, we choose a — 6.042 
in (Q, corresponding to the value S = 1.282 predicted by 
the scaling theory in d = 2 for 6* = —0.282. The plots are 
double-logarithmic plots based on the asymptotic scal- 
ing form ((2Jl with d — 1. It is found, however, that if 
(0) is used as it stands, i.e. L^/^mL{h) is plotted against 
T^i/s^i/*^ it is difficult to get a perfect data collapse at 
small h, for any value of S. This is because the magneti- 
sation per spin in zero field has a binomial distribution, 
and its average over samples, mi(G), is not exactly pro- 
portional to L~^/^ for small L (for L oo, the binomial 
distribution approaches the normal distribution and the 
L^^/^ dependence becomes asymptotically exact). An 
alternative FSS form can, however, be obtained by sim- 
ply replacing the factors L^^^ in (O by l/mi(G). This 
has been done in Fig. 3. In this way, the ordinate is 
identically zero at h — and data collapse at small h is 
greatly improved. 

In the upper plot in Fig. 3, the exact value S — 1.282 
is used, while the lower plot shows the 'best fit by eye' 
obtained with S w 1.37. The curves saturate at large 
h, when all spins are aligned with the field. The scaling 
region extends to, and somewhat beyond, the visible 'el- 
bow' in the data. It is a remarkable that a very good 
fit by eye (lower plot) is obtained for the wrong value, 
d — 1.37, of the scaling exponent, while the exact value, 
(5 — 1.282, is visibly worse (upper plot). 

This paradox, that the best collapse appears to be ob- 
tained with the wrong exponent value, is resolved in a 
similar manner to the h = case. In Fig. 4 we show the 
same data replotted against ln[(L/^)^/^], where ^ is the 
average domain length for a given h, as defined above. 
It is calculated numerically using transfer matrix meth- 



ods and taking systems long enough that the mean do- 
main length converges (which requires longer systems at 
smaller h). 

The collapse is very good, up to and a little beyond the 
elbow, as expected. The disappearance, in the new scal- 
ing variable, of the plateau from Fig. 3 is a consequence 
of the saturation of ^, at the value 2, in this regime, i.e. 
the plateau gets compressed to a single point in these 
variables. As in the zero field case, the failure of the 
data to collapse with the correct scaling exponent, when 
naive FSS is employed, indicates an inappropriate choice 
of scaling variable rather than a failure of FSS itself, i.e. 
the scaling variable L^/'^h^/^ is only useful at very small 
values of h. 

Finally we have studied the 'domain wall' and 'droplet' 
energies directly for the model with bond distribution lpj)l. 
choosing a = 2.546, corresponding to = —0.282. The 
system has L bonds and free boundaries. The domain 
wall energy, E', is given by the magnitude of the weakest 
bond in the system. Averaging over 10^ samples, with 
10 < L < 320, and plotting ln(£') against InL gives the 
exponent 6 = —0.297(2), which is again different from the 
exact value for this model, 9 = —0.282, but the difference 
is much smaller than in the studies described earlier in 
this paper. There is also a small curvature in the data: 
Extracting the exponent from the last two data points, 
L ~ 160 and 320, gives 9 = —0.288, even closer to the 
exact value. 

One can create a droplet excitation by fixing the 
boundary spins in the ground-state configuration and re- 
versing the central spin. The droplet that forms around 
the reversed spin is bounded by the weakest bonds to 
its left and right. Computing the droplet energy numer- 
ically, averaging over 10^ samples for 10 < L < 320, 
and plotting ln(i?dropiet) against InL as before gives 
9 = —0.303(3). Using only the last two data points gives 
9 — —0.290. The discrepancy between the values of 9 
obtained from domain-wall and droplet energies is not 
statistically significant, in contrast to what is observed 
in d = 2 Q . In the latter case, the boundary of a droplet 
forms a closed loop, which may tend to raise its energy 
due to an effective repulsion between different parts of 
the interface, leading to large corrections to scaling p^ . 
In d = 1 the 'interface' consists of two isolated points and 
this effect is absent. 

In summary, we have demonstrated by exact calcu- 
lations in d = 1 that a naive use of FSS, in which 
the asymptotic form of the scaling variable is employed, 
can lead to erroneous estimates of the scaling exponents, 
while FSS itself works rather well. The 'rogue' values of 
9 and S extracted form these studies differ from the true 
values in the same sense as the corresponding exponents 
in d = 2 differ from those obtained using domain-wall 
estimates for 9. We conclude that domain-wall studies 
provide the most reliable determination of the exponent 
9 in spin glasses. 



5 



[1] A. J. Bray and M. A. Moore, J. Phys. C 17, L463 (1984). 
[2] W. L. McMillan, Phys. Rev. B 30, 476 (1984). 
[3] D. A. Huse and I. Morgenstern, Phys. Rev. B 32, 3032 
(1985). 

[4] R. N. Bhatt and A. P. Young, Phys. Rev. B 37, 5606 
(1988). 

[5] N. Kawashima, N. Hatano, and M. Suzuki, J. Phys. A 

25, 4985 (1992). 
[6] H. Rieger, L. Santen, U. Blasum, M. Diehl, M. Junger, 

and G. Rinaldi, J. Phys. A 29, 3939 (1996). 
[7] D. A. Huse and L-F. Ko Phys. Rev. B 56, 14597 (1997). 
[8] M. Pico, F. Ritort, and M. Sales, preprint, 

cond- mat/0210576 (2002). 
[9] K. Binder, Z. Phys. B 48, 319 (1982); W. Kinzel and K. 

Binder, Phys. Rev. Lett. 50, 1509 (1983); Phys. Rev. B 

29, 1300 (1984). 
[10] A. J. Bray and M. A. Moore, in Heidelberg Colloquium 



on Glassy Dynamics, Lecture Notes in Physics 275, 121 
(1986). 

[11] D. S. Fisher and D. A. Huse, Phys. Rev. Lett. 56, 1601 
(1986). 

[12] A. C. Carter, A. J. Bray, and M. A. Moore, Phys. Rev. 

Lett. 88, 077201 (2002); A. Hartmann, A. J. Bray, A. C. 

Carter, M. A. Moore, and A. P. Young, Phys. Rev. B 66, 

224401 (2002). 

[13] N. Kawashima, |cond-mat/9910366| (1999). 

[14] A. A. Middleton, Phys. Rev. B Bs' 060202(R) (2001). 

[15] A. K. Hartmann and M. A. Moore, preprint, 

cond-mat/02 10587 (2002). 
[16] F. Barahona, Phys. Rev. B 49, 12864 (1994). 
[17] W. L. McMillan, Phys. Rev. B 29, 4026 (1984). 
[18] J. D. Reger, K. Binder, and W. Kinzel, Phys. Rev. B 30, 

4028 (1984). 

[19] N. Kawashima and M. Suzuki, J. Phys. A 25, 1055 
(1992). 



